A multi-agent approach to explore the arrival, environment, and adaptation of Homo sapiens in Southeast Asia | Scientific Reports

2021-11-22 08:21:52 By : Ms. Fiona Cai

Thank you for visiting Nature. The browser version you are using has limited support for CSS. For the best experience, we recommend that you use a newer version of the browser (or turn off the compatibility mode in Internet Explorer). At the same time, to ensure continued support, we will display sites without styles and JavaScript.

Scientific Reports Volume 11, Article Number: 21080 (2021) Cite this article

The ability of Pleistocene humans to successfully adapt to different types of tropical forest environments has long been controversial. In order to study the environmental changes in Southeast Asia during the critical period of human species replacement, we analyzed the paleoenvironmental agents of five animal groups in the middle and late to late Pleistocene. The discovery of human teeth has been reported in Duoi U'Oi (70-60 ka) in Vietnam and Nam Lot (86-72 ka) in Laos. However, the use of paleoproteomics has allowed us to abandon the latter, and so far, no human remains over 70 ka have been recorded in this area. Our research results show that tropical rain forests are highly sensitive to climate change during that period, and canopy forests fluctuate greatly. Locally, large fauna is resilient to these fluctuations, until the cooling period of marine isotopic stage 4 (MIS 4; 74-59 ka) changes the entire biological community. Then, under strong selection pressure, populations with new phenotypic characteristics appeared, while some other species disappeared. We believe that this climate-driven transition provides humans with new foraging opportunities in the new rainforest environment, and is likely to be a key factor in the settlement and spread of our species during MIS 4 in Southeast Asia.

Although the earliest forms of humans occupy different C3-C4 environmental niches in Africa1, this genus is generally considered to be primarily adapted to open environments2,3. In Asia, early Homo erectus may have lived in areas without forests along the banks of northern China and the Java River valley, and the niche division between ancient humans and other large primates living in densely forested habitats has been proposed4,5, 6,7,8.

During the late Pleistocene, with the extinction of the last Homo erectus in Indonesia8, the last Denisovans may have existed in several parts of the continent9, and eventually all ancient groups were replaced by Homo sapiens upon arrival. 10. In On the continental scale, it has been suggested that the transition from open habitats (mixed savanna and woodland) to rainforest habitats during the transition period between the Middle Pleistocene and Late Pleistocene triggered the decline of ancient humans and was unable to adapt to these new environments11.

Therefore, determining the paleoenvironmental background faced by different human species in Southeast Asia (SE) Asia may trigger debates about the uniqueness of our species. However, in the Pleistocene tropical Indochina subregion, rare dental remains temporarily attributed to ancient humans are often reinterpreted as the remains of great apes (mainly orangutans of the genus Orangutan)5,12,13,14, so it is difficult to associate ancient humans with The records are firmly linked to past vegetation in many cases.

The propagation route of H. sapiens to southern China may have crossed Indochina15, but the time of this event, its course—maybe one or several waves (ka)16,17 about 100,000 years ago—and The problem of how H. sapiens adapts to the rainforest environment remains unresolved. Of course, the lack of detailed chronology of several Southeast Asian sites has caused our understanding of this period to become blurred. So far, the earliest indisputable archaeological evidence of human adaptation to Asian tropical rain forests is actually the most recent, dating back to Lida Ajel, Sumatra, about 73-63 ka.

Here, we try to solve some of the key issues by analyzing five mammal groups from Vietnam and Laos to understand the evolution of our species. Their age ranges belong to different marine isotopic stages (MIS): Coc Muoi (148 – 117 ka, MIS 6–5), Tam Hang South (94–60 ka, MIS 5–4), Nam Lot (86–72 ka, MIS 5), Duoi U'Oi (70–60 ka, MIS 4) And Tam Hay Marklot (38.4–13.5 ka, MIS 3–2) 19, 20, 21 (Figure 1). From a paleoecological point of view, the crown size and stable isotope measurements of taxa identified from these faunas are agents of environmental reconstruction7,11,20,22,23,24 and organisms occupied by humans The main source of information for the group. In the area studied, the earliest records of Homo sapiens appeared in Tam Pà Ling (~ 70–46 ka25,26) and Duoi U'Oi (70–60) with several human skeletal remains of approximately 70 ka (ka22). However, the older hypothetical ancient human specimens associated with the Nam Lot combination (86-72 ka22) opened up the possibility of earlier arrival27. Therefore, we use Paleoproteomics 28, 29, 30 to solve the specific allocation of this specimen based on its enamel proteome, with the purpose of better placing the arrival of modern humans in the local context.

(A) Location of locations in northern Laos (Tam Hang South, Nam Lot and Tam Hay Marklot) and northern Vietnam (Coc Muoi and Duoi U'Oi). (B) The δ18O records of the Chinese caves in Sanbao, Dongee and Hulu show millennium-scale climate changes (modified after 87) related to the 224 ka East Asian summer monsoon intensity changes. The decrease in the δ18O value (‰, VPDB, left vertical axis) corresponds to an increase in precipitation, which is a quantitative effect 77. The right vertical axis corresponds to the summer sunshine in the northern hemisphere (65°N, W m-2). The age range of the fauna is placed below the δ18O record curve, from the oldest (right) to the youngest (left).

From the middle to late Pleistocene, the local fauna consisted of most modern taxa and some ancient taxa 23, 31, 32, 33, 34, 35. In general, these faunas are similar to those on other continents at the time, because they were dominated by large herbivores (> 1000 kg to 5000 kg36), including elephants, saber-toothed tigers, giant tapirs, and several rhinoceros and large bovids. Animal 37. This similarity also extends to the local associations of species with different ecology and the differences in their spatial distribution from today (for example, orangutans and giant pandas or tigers and hyenas). Most paleontologists now believe that this unexpected association of species named dissimilar fauna (38,39,40,41,42, but see Reference 43) is the result of species’ different responses to environmental changes. That is in the "individualistic way"-according to their life history 44, 45, 46, 47.

Using classical animal archeology methods, the paleoecological reconstruction analysis of the Southeast Asian mammalian community failed to find the main functional changes of the mammalian community 37, 48, 49, and paleontological records seem to be somewhat consistent. In fact, in this tropical region, this species shows a wide ecological range, including latitude and altitude. Nevertheless, the question remains how species survived climate change during the Pleistocene and adapted to non-similar climates, exhibiting different sets of environmental variables (vegetation structure, sunshine rate, seasonality, rainfall, etc.). Today's 50,51. In addition, the evolution of various lineages at the subspecies level is usually unknown, because the movement of populations can rarely be tracked in the fossil record. Another limitation of the study of species dynamics in tropical Asia is that there is no preserved DNA in the fossil remains, which hinders the reconstruction of the genetic history of mammalian lineages. In Eurasia, most of the molecular analysis of ancient DNA (aDNA) focused on the Bering fauna, emphasizing that the Late Pleistocene was a dynamic period of cold-adapted mammals affected by climate change 52,53,54,55,56 ,57,58,59, 60, 61, 62, 63, 64. These studies prove that various processes (ie, population shrinkage, local extinction, migration, replacement of new populations, or interspecies competition) lead to the success of new and better-adapted populations over time 47,65,66,67. At the same time, the impact of a colder climate on people who adapt to warmth remains largely unknown.

However, the use of morphometry, especially the isotope replacement of teeth, can solve various environmental problems, especially the assessment of the impact of large-scale climate oscillations on the regional-scale rainforest ecosystem and its impact on mammalian communities and related humans, such as Yue Lai The more research 20, 22, 24, 68, 69, 70, 71, 72, 73.

Our data set contains six mammalian orders, namely Artiodactyla, Proboscis, Carnivora, Primates and Large Rodents (Methods, Supplementary Materials and Methods, Supplementary Figure S1-S3, Supplementary Table S3-) Hundreds of isolated teeth of mammals S4). All stations are located in a narrow latitude zone between 23° and 20° passing through Laos and northern Vietnam (Figure 1A). Due to the abiotic parameters related to the latitude distribution, namely climatic effects (temperature, distance from the coast, seasonality of rainfall, rainfall), the location of these locations minimizes changes in the size of the species. However, these five locations are located at different altitudes, from lowland locations in the alluvial plain (Duoi U'Oi) to Zhongshan locations (Nam Lot and Tam Hang South).

First, we used new data from Coc Muoi, Duoi U'Oi, and Tam Hang South to compare carbon (δ13Capatite) and oxygen (δ18O) isotopic measurements from tooth enamel. These data came from 335 taxa belonging to a wide range of taxa. Specimens, as well as those published in Nam Lot and Tam Hay Marklot20,22. We estimated the δ13C carbon source value in animal diets to specifically analyze the changes in the ratio of C3 plants (trees, shrubs, shrubs, and grasses) to C4 plants (grass, sedges) during the study period74. The δ13C of biological apatite allows the reconstruction of the paleoenvironment based on these isotope-different carbon sources. We also use δ13C carbon sources to separate C3 canopy forests from other C3 forest habitats to reveal local fluctuations associated with climate change in tropical rain forests75,76. The δ18O value is used to provide additional paleoecological information related to changes in abiotic conditions (latitude, climate, temperature, moisture content, precipitation, and isotopic composition) 77, 78, 79, 80, 81, 82, 83, 84 ("Methods ").

In addition, we used raw morphological measurement data-the crown area of ​​213 specimens belonging to five taxa of herbivores and omnivores, sambar (Rusa unicolor), muntjac (Muntiacus sp.), mountain lion (Capricornis sumatraensis) , Wild boar (Sus scrofa) and macaque (Macaca sp.)-to detect significant phenotypic changes in their pedigree over time. Combining proxy data based on stable isotope data with these morphometric data allows us to determine which climate changes have the greatest impact on mammalian communities associated with rainforest dynamics.

Our fourth agent considers the type of digestive physiology and uses the ratio of ruminants to hindgut fermenting herbivores by weight category as an indicator of open landscape expansion (mainly through the emergence of exclusive grazing groups), so rainforest .

Finally, we discussed how the climate change that occurred during the late Pleistocene affected the local adaptation of the first batch of Homo sapiens, as well as the wider adaptation in Southeast Asia. To this end, we used available climate records, such as pollen data 85 and 86 and Chinese cave δ18O data 87 from cave animals, as other relevant information sources.

The MS/MS spectrum clearly assigns Nam Lot incisor (NL 433) to the genus Orangutan (orangutan), which has no unique and highly reliable match with Homo29,30. For those positions where we have the proteomics data of Nam Lot specimens, in our reference sequence, there is no sequence difference between Pongo abelii and P. pygmaeus. Therefore, we assigned the specimen to the Pongo genus without further specifying the species (additional methods and results).

The δ13Csource and δ18Oapatite values ​​of specimens belonging to all taxa are shown in Figure 1 and Figure 2. 2 and 3 and supplementary tables S5-S7. The values ​​of all samples and reference standards (δ13Capatite, δ13Ccarbon source and δ18Oapatite) are shown in Supplementary Appendix S1-S2. For the new data measured here, the δ13C carbon source values ​​of Coc Muoi, Duoi U'Oi and Tam Hang South range from-33.8 to-18.1 ‰ (average δ13C carbon source =-28.0 ± 2.4 ‰ (1 σ), n = 84 ), − 34.3 to − 15.1 ‰ (average δ13C carbon source = − 28.4 ± 3.0 ‰ (1 σ), n = 60) and − 30.0 to − 12.0 ‰ (average δ13C carbon source = − 25.0 ± 3.6), (‰ 1 σ) n = 62), respectively (Figure 2).

The histogram distribution of the relative frequency (%) of the δ13C carbon source values ​​of all taxa in the five Southeast Asian faunas, arranged in chronological order from the oldest (left) to the youngest (right). Each bin represents an interval of 1‰. The shaded area represents the δ13C carbon source value related to the closed forest canopy (δ13C carbon source<-27.2‰); intermediate rainforest and woodland biological communities (δ13C carbon source> − 27.2 ‰ and <− 21.3 ‰; and savannah environment (δ13C carbon source) > − 15.3 ‰). The white area (δ13C carbon source> − 21.3 ‰ and <− 15.3 ‰) is the value derived from the comprehensive consumption of C3 and C4 resources, and does not correspond to any specific ecological environment. The red dotted line represents the δ13C carbon source of each site Value mean.

The distribution of δ18O values ​​of all taxa in the five Southeast Asian faunas, in chronological order, from oldest (left) to youngest (right): Coc Muoi (), Tam Hang South (), Duoi U'Oi () and the former Nam Lot () and Tam Hay Marklot () released data. The outline of the violin chart represents the nuclear probability density, where the width shows the proportion of the data found there. The box in the box and whisker diagram in the violin chart represents the 25th to 75th percentile, and the median is a thick horizontal line.

The new δ18Oapatite values ​​obtained at the three sites range from -9.8 to -3.9 ‰ (average δ18Oapatite =-6.1 ± 1.2 ‰ (1 σ), n = 84), -9.7 to-2.5 ​​‰ (average δ18Oapatite = -1.6. ‰ ( 1 σ), n = 60) and −9.4 to −2.8 ‰ (average δ18Oapatite = − 6.7 ± 1.4 ‰ (1 σ), n = 62), respectively Coc Muoi, Duoi U'Oi and Tam Hang South (Figure 3 ).

The statistically significant difference between sites from novels (Coc Muoi, Tam Hang South and Duoi U'Oi) and published data (Nam Lot and Tam Hay Marklot) is the source determined by one-way analysis of Kruskal-Wallis δ13Ccarbon variance (H(4) = 83.3, p value <2.2e-16) and δ18Oapatite (H(4) = 25.5, p value = 4.019e-05). Afterwards, Dunn's pairwise comparison of the test showed that the δ13C carbon source values ​​of the Vietnam site and the Laos site are different. Tam Hay Marklot and Nam Lot also seem to be very different. Finally, the δ18Oapatite value from Tam Hang South was determined to be significantly different from the values ​​of all other sites except Duoi U'Oi, and Duoi U'Oi itself was extracted to be significantly different from Nam Lot and Tam Hay Marklot ( Supplementary Table S8, S9).

Broadly speaking, the range and median fluctuation of the δ18O value (Figure 3) are consistent with the δ13C carbon source value (Figure 2). However, the distribution of δ13C carbon source values ​​highlights that during the study period, the C3 forest environment (canopy forest, middle rain forest and woodland) still dominates. In addition, when we look at the percentage of samples based on the distribution of δ13C carbon source values ​​related to different biomes in Table 1, the data shows that environmental conditions have changed significantly through the Coc Muoi-Tam Hay Marklot time series. Therefore, tropical forests are obviously very sensitive to climate change. Our research results specifically illustrate the dynamics of canopy forests (δ13C carbon source <− 27.2 ‰), and show their potential to shrink across time and space: Coc Muoi (65.4%), Tam Hang South (27.4%), Nam Lot ( 42.1%) ), Duoi U'Oi (73.3%) and Tam Hay Marklot (26.3%).

The sequence of animal groups by body weight and digestion strategy is shown in Figure 4 and Supplementary Table S10. Among the three oldest animal groups, the hindgut fermentation group, that is, the non-ruminant group, includes seven large herbivores (> 350 kg) and large herbivores (> 1000 kg), belonging to the following genus: giant tapir , Tapir, saber-toothed tiger, elephant, rhinoceros and Dicerorhinus (compared to only one ruminant Bos species) dominate the biomass. Duoi U'Oi’s “ruminant and non-ruminant taxa” has a ratio of 4:7, indicating that the composition of giant herbivores has changed, without giant tapirs (> 350 kg) and saber-toothed elephants (> 1000 kg) . However, hindgut fermented herbivores still dominate, because the increase in ruminants cannot compensate for the loss of diversity in the huge ancient taxa. Tam Hay Marklot marks a shift, representing small and medium ruminants (18 to 350 kg) (including Rucervus eldii, Axis porcinus and Naemorhedus caudatus) becoming dominant (ruminant to non-ruminant group ratio is 8:6). This trend has clearly continued to the present, as seen by the increase in grazing species in the current fauna at these latitudes (ratio 9:4) (Figure 4).

The number of species in the five animal groups divided by body weight category and digestion strategy, from the oldest (left) to the youngest (right) in chronological order. This ratio refers to the number of ruminant and non-ruminant taxa. For a list of taxa in each weight category, see Supplementary Table S10.

The size range of the crown area of ​​p3 (Rusa unicolor and Sus scrofa) and m3 (Macaca sp., Muntiacus sp. and Capricornis sumatraensis) varies by location (Figure 5A). However, statistical analysis is limited by the imbalance of sample sizes in certain locations. Only R. unicolor (n = 65) and S. scrofa (n = 61) use the Kruskal-Wallis test (H(4) = 21.09, p-value = 0.0003 and H(4) = 14.25, respectively, p-value = 0.007). Afterwards, Dunn's test paired comparison extracted from Coc Muoi R. unicolor and Nam Lot (p value <0.005) and Duoi U'Oi (p value <0.05) were significantly different; and R. unicolor's Nam Lot and Tam Difference in Hay Marklot (p <0.05). S. scrofa samples also showed significant differences between populations (p value <0.05): Coc Muoi and Nam Lot; Tam Hang South vs. Duoi U'Oi; Nam Lot and Duoi U'Oi and Tam Hay Marklot (Supplementary Table) S11 and S12).

The distribution of crown area size (A) and δ13C carbon source value (B) in the five taxa ) To the youngest (bottom) in chronological order. CM Coc Muoi, THS Tam Hang South, NL Nam Lot, DU Duoi U'Oi, THM Tam Hay Marklot. See Supplementary Table S13 for sample size.

In general, there is a cumulative effect from Coc Muoi to Nam Lot, that is, the population follows directed evolution, towards larger (R. unicolor, S. scrofa and Macaca sp.) or smaller (Muntiacus sp. and C. sumatraensis) The surface of the crown area is based on taxa. In the entire fauna sequence, Duoi U'Oi marks this shift in directed evolution. This shift is particularly pronounced in S. scrofa, but the five taxa studied seem to be affected by this reversal of pedigree dimensions (Figure 5A). This reversal is used here as a signal that, in the face of high selection pressure, due to population adaptation or replacement (through extinction or assimilation), populations with new phenotypic characteristics have emerged. Therefore, in the Coc Muoi-Tam Hay Marklot time series, Duoi U'Oi seems to indicate a certain population turnover.

Judging from the available records discussed here, there was no obvious human presence in this area before ~70 ka. However, as shown by our paleoproteomic analysis of Nam Lot incisors (86-72 ka), the ability to obtain protein sequence information from tropical regions and distinguish between Pongo and Homo opens up the possibility of directly solving this problem in the future in Southeast Asia The early Homo sapiens that existed.

The relative similarity of δ13C carbon sources and δ18O values ​​between Coc Muoi (148–117 ka) and Duoi U'Oi (70–60 ka) indicates that climatic conditions have induced an ecosystem dominated by C3 in two different periods. As shown in the curve of Sambo/Calabash δ18O Chinese Cave Record 87 in Figure 1B, considering the age range of the fauna, the advantages of these forest ecosystems may be related to the two high-level declines in monsoon intensity during MIS 6 Coc Muoi (MIS 6.288 ,89) and Duoi U'Oi's MIS 4 period (Figure 2 in87). In Coc Muoi and Duoi U'Oi, the enclosed rainforest contains most of the mammalian biomass, mainly composed of browsers weighing up to 5000 kg (Figure 4 and Supplementary Table S10). However, upon closer inspection, the two locations revealed significant differences in species that rely on the canopy forest as a diet (Figure 6). First, Duoi U'Oi marks a decline in the diversity of large herbivores, because there are no two ancient taxa: the giant tapir Megatapirus augustus and the proboscis saber-toothed tiger. Both sites are located in the same vegetation area, with an altitude of <400 m (asl), and other sources of variation are reduced, supporting the hypothesis that climate has a major impact on mammalian communities.

Compare the δ13C value ranges of selected taxa in the five animal groups in chronological order from the oldest (left) to the youngest (right). The shaded area represents the δ13C carbon source value related to the closed forest (δ13C carbon source <− 27.2 ‰), intermediate rain forest and woodland biome (δ13C carbon source> − 27.2 ‰ and <− 21.3 ‰); and the environment similar to the tropical grassland ( δ13C carbon source> − 15.3 ‰). The white area (δ13C carbon source>-21.3 ‰ and <-15.3 ‰) is the value of the comprehensive consumption of C3 and C4 resources, and does not correspond to any specific ecological environment.

Second, the results of the δ13C carbon source value show that in Duoi U'Oi, environmental changes have led to the redistribution of niches and interactions between new species. The difference from the old Tam Hang South and Nam Lot fauna is significant, and the MIS 4 Duoi U'Oi fauna seems to be significantly different from the MIS 5 fauna (Figure 6 and Supplementary Table S8). Especially the sambar R. unicolor, tapir and large bovids that forage in this new Duoi U'Oi ecosystem dominated by C3. In addition, the Duoi U'Oi δ18O value is lower than the value prevailing in Nam Lot globally, which indicates that colder and/or wetter conditions are present (Figure 3).

Although the forest environment dominates during this period, some animal groups also include species that rely on mixed C3-C4 and/or C4 resources for their diets (Figures 2 and 6). Especially in the case of Tam Hang South and Nam Lot fauna, the δ13C carbon source value shows the existence of a more open environment and the reduction of forest habitat: Tam is 11.2% (C3–C4) and 3.2% (C4) Hang South and Nam Lot They are 5.2% (C3–C4) and 3.5% (C4) respectively (Table 1). These data are related to seasonal increases24. Zheng Helei's data (from near the Leizhou Peninsula in southern China, latitude 21°-20°, elevation <260 m85) shows that the mountain slopes are covered by monsoons with evergreen forests and dense shrubs. Tam Hang South and Nam Lot also show different sets of environmental variables, and their δ18O values ​​are significantly different (but not due to altitude changes, because they are about 150 m apart and located at the same altitude ~ 1120 m), indicating that despite Lack of age accuracy, how dynamic is the interglacial MIS 5?

In the second half of the Late Pleistocene, there were no records of mammalian fauna in the study area. The carbon and oxygen isotopic records of land snails indicate that between ~70 and ~32 ka, the landscape is still highly forested and humid. Camaena massiei is located near Laos The site of Tam Pà Ling 90. The significant changes in the more open landscape associated with increased drought due to low sea levels and the expansion of the continental land surface were later illustrated by our Tam Hay Marklot (38.4-13.5 ka20) δ13C carbon source value (Figure 6). Further west, the Tham Lod Rock Sanctuary in Thailand (33–11.5 ka70,92) shows a similar situation: mixed breeders, sambar and large bovines seem to have migrated to grassland areas and herbivores and medium-sized deer. (Axis, Rucervus) and calves (Naemorhedus) occupy new spaces, most likely from southern latitudes. Over time, these subpopulations have adapted to living in open landscapes93 (Supplementary Figures S5 and S6) , Supplementary Appendix S4). Boh Dambang (25-18 ka, Cambodia 94), 91 further south in the center of the hypothetical savanna corridor and the site of Tam Hay Marklot show that despite this, the peninsula is still dominated by herbivores (Muntiacus, Rhinoceros, and tapir), even though During the last glacial maximum (LGM; 26.5-19 ka) the drought height near 95. In fact, Tam Hay Marklot's δ13C carbon source value of 26.3% and Boh Dambang's 9.6% indicate this environment (Supplementary Table S14), and the results are from a fairly comparable sampling classification group at the two locations (Supplementary Table S7). The shrinkage of forest habitats has significantly changed the carrying capacity of the ecosystem, forcing large carnivores to look for new hunting opportunities among species that inhabit the savanna, including tigers (Panthera tigris), leopards (P. pardus) or hyenas (Crocuta crocuta) ) (Figure 6 and Supplementary Figure S5).

The current analysis combines the canopy area size (Figure 5A) with the δ13C carbon source data (Figure 5B), revealing that the main turnover and decline of the population in the Duoi U'Oi mammal community is related to the regression of the C3 dominant landscape. Isotope value. In Duoi U'Oi, populations with new phenotypic characteristics have emerged compared to the old locations, and they have adapted better to this new environment. This series of evidence can be interpreted as an adaptive response to major selection pressures. In fact, species of different weights range from ~ 15 kg (muntjac) to ~ 220 kg (sambar), different dietary strategies (omnivores and herbivores), and adapt to different ecological niches (ground-dwelling ungulates and trees). Perching monkeys) ), experienced similar evolutionary trends, indicating that the entire ecosystem has been affected. As shown in Figure 1B, the curve recorded by the Sambo/Gourd δ18O cave shows that this period witnessed short-term climate changes, and the monsoon intensity dropped sharply at the beginning of MIS 4 (events of comparable magnitude occurred during the LGM period). Such climatic events may lead to strong selection pressure, which triggers new adaptations and population migration. The suddenness of this climate change lasting about <300 years and the average temperature drop of about 5-6 °C85 are the most likely explanations for the population turnover and the disappearance of the last ancient species (sabertooth and sabertooth). Giant tapir). So far, from the animal records studied, there is no evidence that these species are related to Homo sapiens. Increasing evidence based on systematic geographic analysis (aDNA) indicates that population replacement within a species may be a rapid process leading to major extinction and recolonization events due to interactions between populations or due to the impact of sudden climate change 53, 54.

In southern China, the late Middle Pleistocene to early late Pleistocene series of Yugong, Quzhai, and Baxian show the same trend, as shown by the δ13C value patterns of R. unicolor and S. scrofa24,35,73 (Supplementary Figure S7 and S8, supplementary appendix S3). Therefore, the cooling event may lead to the expansion of south China canopy forest dwelling populations to low latitudes.

From a paleoecological perspective, our findings confirm that tropical rainforest ecosystems prevailed in the late Mid-Pleistocene at these latitudes (Supplementary Figure S7). They are consistent with the environmental reconstruction done by Louys and Roberts. Our findings also question the relationship between the changes in the tropical rain forest environment and the main turnover of early humans (Homo erectus, Denisovan) in Southeast Asia before the arrival of Homo sapiens. In fact, assuming that ancient humans could not adapt to tropical rainforest habitats5,11,17, environments like Coc Muoi might become obstacles to their local settlement.

As our research shows, environmental changes not only affect the fauna, but may also affect Homo sapiens. In fact, in the extensive record of forest persistence in the late Pleistocene, we did notice important changes in the area when it was first known to occur. The two severely worn teeth of Duoi U'Oi (70–60 ka22) (thus being carefully classified as ancient or human) recorded the existence of humans, while Homo sapiens appeared in Tam Pà approximately 200 kilometers away Ling cave ruins time (~ 70–46 ka25,26). Judging from the existing records, it is clear that humans settled in the area around the time of MIS 4. 70-50 ka97, must deal with densely forested environment 90. However, the relative cooling period of MIS 4 may have led to a profound transformation of forest composition and structure, as evidenced by Zheng Helei's 85 pollen record analysis (altitude <260 m) at the same latitude and further east. In fact, due to the decline of the mountain vegetation zone, the monsoon evergreen forest is transformed into a mountain forest, and temperate plants have increased, and conifers, previously rare trees and ferns (20%) have increased significantly. According to the pollen records of the caves in southern China at a comparable altitude (< 212 m), Li et al.86 emphasized similar changes in the late Pleistocene subtropical mixed coniferous forest.

Although the palynological analysis of Duoi U'Oi is based on few elements, it also shows that the frequency of fern spores (25%) and non-arboreal groups (31%) is relatively high, but the representativeness of mangrove pollen grains is relatively high. Low (2%) 22. In the absence of archaeological evidence (stone tools or organic industries, bones with traces of slaughter, etc.), it is still challenging to assess many aspects of human foraging behavior in this region of Southeast Asia. However, different types of tropical rain forests provide edible plants, fruits, seeds, nuts, and honey98, while rivers provide predictable resources of shells, fish, molluscs, and algae. Despite the loss of the ancient large herbivores, the Duoi U'Oi biota supports a wide range of games. The presence of humans on Duoi U'Oi is recorded not only through two isolated teeth related to the animal group, but also through indirect evidence. In fact, the mortality characteristics of sambar deer at this location clearly show human characteristics19. In Duoi U'Oi, humans can hunt down prey deep in the rainforest and selectively hunt adult individuals. It is possible that during MIS 4, new vegetation that changed shrubs, ferns, and herbaceous strata made the forest easier to enter and navigate, and provided new hunting opportunities for foragers. Mankind has clearly successfully adapted to the general habitat of the closed canopy forest 18,99,100,101, probably long before the formation of specialized foraging behavior, especially in the hunted prey, arboreal and semi-arboreal species and terrestrial creatures. The proportion of species is larger ~45 ka99, 102,103,104.

Due to the climate dynamics of the late Pleistocene, our research results support the successful spread of our species to the area during MIS 497, 105, 106, and 107. The combination of ecological and behavioral factors seems to help early Homo sapiens successfully cope with the challenges of the rainforest environment: the renewal of vegetation allowed humans to occupy a new niche at the beginning of MIS 4108, and our ability to adapt species effectively to this environment109 . In view of the recent reassessment of human spread in the Far East. From 65 ka to 45 ka, despite the great controversy 97,110,111, our species may enter southern China through similar types of rainforests. The similarity of isotope data between Duoi U'Oi and Baxian24,35 also shows (Supplementary Figure S7) ). In the middle and late Pleistocene, there was no comparable isotope data available on the Southeast Asian continent, but it remains to be proved that Homo sapiens expanded its habitat during this period when it moved to the Southeast Asian islands. However, we know that with this climate-driven vegetation replacement, modern humans arrived in Sumatra in 73-63 ka, where they occupied the dense evergreen rainforest ecosystem18. Starting at 46 ka, in the Niah Cave in Borneo, modern humans have effectively used the tropical environment, used sophisticated hunting techniques and were able to process poisonous plants for food103. Interestingly, climate change during MIS 4 is considered to be the main driving factor for human migration from Africa to Eurasia. The reason is that the colder weather and lower sea level 112,113,114 make it possible for the population to reach northern Australia before 65 ka115,116. .

Studies on the evolution of ecosystems in Southeast Asia in the past have been limited by the scarcity of paleontological records, atypical remains protection, discontinuity, and the accuracy of the available age of the fauna, which limit our understanding of the warm-adapted species response to climate change. Nonetheless, our multi-agent method combines morphometric and isotope data as well as digestive physiology types, revealing that the cooling event of MIS 4 may have a profound impact on the overall biological community of the area when humans arrive, resulting in a relationship with dense forests. Landscape-related mammal populations. Our research results indicate that the new composition and structure of tropical rain forests are likely to be key factors in promoting the rapid spread of Homo sapiens to Southeast Asia.

In addition, our analysis emphasizes the importance of paleoproteomics in elucidating the taxonomic distribution of remains (paleohumans and equines). In fact, given the scarcity of evidence of early Homo sapiens in the region, this method is essential for further understanding of the arrival and spread of our species on the continent.

These stations are located at different altitudes: 113 m (asl) above sea level (Duoi U'Oi, Vietnam), 361 m asl (Coc Muoi, Vietnam), 1120 m asl (Tam Hang South/Nam Lot, Laos), and 809 m asl (Tam Hay Marklot, Laos). They are located 120 kilometers (Duoi U'Oi), 170 kilometers (Coc Muoi) and 270 kilometers (Lao ruins) from the coast bordering the Gulf of Tonkin (Supplementary Figure S1). Today, the region is characterized by a humid subtropical climate. According to the Köppen-Geiger climate classification, summers are hot and humid, and winters are cold and mild. Local climate data, annual average temperature and annual average rainfall are: 23.7 °C and 1735 mm (Hoà Binh Province, Duoi U'Oi, Vietnam); 21.9 °C and 1349 mm (Lang Son Province, Coc Muoi, Vietnam) ; 19.8 °C and 1331 mm (Houaphan Province, Xamneua, Nam Lot, Tam Hang South and Tam Hay Marklot, the highlands of Laos) (https://en.climate-data.org/).

These sites were excavated between 2003 and 2015. Descriptions of the site and paleontological content can be found elsewhere in the publication 19, 20, 21, 22 and condensed versions of supplementary data. Due to the burial and geological processes of karst deposits, all combinations show similar remains preserved. They are mainly composed of isolated teeth of mammals, and due to preservation deviation, only species composition combinations> ~ 5 g: Coc Muoi (NISP = 1323), Tam Hang South (NISP = 673), Nam Lot I (385 teeth) and Nam Lot II (5 teeth) (NISP = 390), Duoi U'Oi (NISP = 871) and Tam Hay Marklot (NISP = 1364) (Supplementary Tables S3 and S4). Soil analysis of these sites shows that porcupines are the main accumulator of large mammal bones (most of their teeth are bitten) before they are buried in sediments. All the remains deposited in these locations were flowed through the entire karst network with water that caused the loss of the smallest elements19. According to field observations, Brain117 shows that the remains collected by porcupines can well represent a place or the number of carcasses around a place, so it is a good representative of species richness. Therefore, these combinations consist of similar taxa. Almost all groups can sample for stable carbon and oxygen isotope analysis (except Proboscidea from Tam Hang South and Tam Hay Marklot), allowing comparison of communities over time (Supplementary Table S7).

These five animal groups do not represent continuous records, because there are two main gaps, namely Coc Muoi (148–117 ka) and Tam Hang (92–60 ka) and Duoi U'Oi (70–60 ka) and Marklot (38.4) –13.5 ka) (Figure 1 and Supplementary Figures S2 and S3). No artifacts or other objects (charcoal, decorations, traps, etc.) related to the remains of the animals were found.

In this analysis, we use the size of the crown area of ​​the tooth (maximum length x maximum width) as an indicator of ecological changes in five mammalian lineages (Muntiacus sp., Capricornis sumatraensis, Rusa unicolor, Sus scrofa, and Macaca sp.). We chose these taxa defined at the species or genus level because they are common to all five animal groups and have a sufficient number of specimen records (Supplementary Tables S5 and S6 and Annex S6). However, this analysis is limited by many biases, including the differential representation of tooth types within a given taxon. This is why we choose left and right m3s among Muntiacus sp., C. sumatraensis and Macaca sp. Regarding the other two taxa, the number of several tooth types is significant, p3, p4, and m3 in S. scrofa, and p3 and p4 in R. unicolor, but they have different ranges of variation. In R. unicolor, unlike p3s, p4s has a similar range of variation between sites. In S. scrofa, m3s shows greater variability, and there is a large overlap between sites. Considering the purpose of our research, using the crown area size as a signal for the new phenotype p3 is the most useful tooth type in the two taxa. This is most likely due to selective adaptation pressure on the skull. For example, in suids, the populations differ in the length of the skull, upper jaw, and row of teeth118. We did not try to estimate the weight of an individual.

The stable carbon isotope of bioapatite (δ13Capatite) reflects the relative proportion of carbon in the consumer's diet, which comes from the main source of the food web, that is, plants that use the C3 or C4 photosynthetic pathway74. In tropical and subtropical regions, wetter forests and woodland habitats are associated with C3 plants that exhibit low δ13C values, while drier and more open environments are characterized by C4 plants with high δ13C values76,119,120. In addition, the lowest δ13C value reflects the dense forest conditions caused by the "canopy effect", so it can be used to distinguish the C3 forest environment. Finally, using the measured δ13Capatite value and the enrichment factor adapted to body weight (supplementary materials and methods), we estimated the initial δ13C value of the carbon source in the animal's diet, which is labeled "δ13C carbon source" here.

Stable oxygen isotopes (δ18O value) are used to provide additional paleoecological information. The main source of δ18O changes in enamel is the oxygen isotope composition of chemically bound water (ie water found in plants) in drinking water and diet 78,79,80,81,82,83. This kind of water itself is controlled by various environmental and geographical conditions such as latitude, climate, temperature, moisture content, precipitation, precipitation isotope composition, etc. (In low latitude areas, the change of δ18O precipitation is mainly affected by precipitation, that is, the effect) 77,78 ,84.

Assigned to orders of mammals from Coc Muoi (n = 84), Tam Hang South (n = 62) and Duoi U'Oi (n = 60) (Artiodactyla, Proboscis, Carnivora, Primates and Rodents) ) The fossil teeth were sampled and analyzed for this study (Supplementary Tables S5 and S7). The enamel is first mechanically cleaned using a hand-held dental drill equipped with a diamond grinding head. Use diamond tip cutting wheels or diamond tip burrs, and then collect samples (powder or fragments) along the entire crown height of each sample. When sampling enamel fragments instead of powder, use an agate mortar and pestle to crush the complete enamel fragments. The enamel powder sample was soaked in 1 ml CH3COOH (0.1) M for 4 hours at room temperature, then rinsed several times in distilled water, and finally dried overnight at 65 °C. Using the carbonate phase of the enamel, in the "Service de Spectrométrie de Masse Isotopique du Muséum (SSMIM)" in Paris, a stable carbon and oxygen isotope ratio was measured using a Thermo Scientific Delta V Advantage isotope mass spectrometer and Thermo Scientific Kiel. IV Carbon Salt device chemicals. Isotopic abundance is expressed in δ (delta) symbol, expressed as deviation per mil (‰), where: δ13C = (13C/12Csample/13C/12Cstandard − 1) × 1000 and δ18O = (18O/16Osample/18O/16Osta​ ​Ndard − 1) × 1000.

In this analysis, we used the δ13C limit corresponding to a wide range of pre-industrial environments:-27.2 ‰ and-21.3 ‰ as the upper limit of δ13C for closed canopy forest 121 and intermediate rainforest and woodland biological communities, which are 74, 122 and-respectively 15.3 ‰ is 74 as the lower limit of δ13C for the C4 savanna environment.

Kruskal-Wallis one-way analysis of variance was performed on the data sets of novels (Coc Muoi, Tam Hang South and Duoi U'Oi) and published data (Nam Lot and Tam Hay Marklot) to determine the statistical difference and location of δ13C carbon sources. The value of δ18Oapatite between points. To this end, samples come from Coc Muoi (n = 84), Tam Hang South (n = 62), Duoi U'Oi (n = 60), and Nam Lot (n = 5722) and Tam Hay’s published sites use Marklot (n = 7220). The size of the crown area between the stations was also studied for R. unicolor (n = 65) and S. scrofa (n = 61). All analyses chose the Kruskal-Wallis test instead of the parametric analysis of variance, because a preliminary test was performed to check for normally distributed data and equal variances, which indicates that a nonparametric test will be used. All statistical analysis was performed using the free program R software (R Core Team, 2018).

Sponheimer, M. Isotopic evidence of early human diet. The process of the national team Akkad. science. United States 110, 10513–10518 (2013).

CAS PubMed Central Article ADS PubMed Google Scholar 

Fleagle, JG, etc. (eds) Out of Africa I: The first ancient human colonization of Eurasia. Vertebrate Paleontology and Paleoanthropology (Springer, 2010).

Norton, CJ & Braun, DR (eds) Paleoanthropology of Asia: From Africa to China and other regions. Vertebrate Paleontology and Paleoanthropology (Springer, 2010).

Betis, EA III. Wait. Out of Africa: Early Pleistocene environment inhabited by Homo erectus in Sanjilan, Java. J. Humph. evolution. 56, 11-24 (2009).

Ciochon, RL Separation of ancient humans from the sabertooth-giant panda fauna: new perspectives on ancient Asian humans. In "Out of Africa I: The First Ancient Human Colony in Eurasia" (eds Fleagle, JG et al.) 111–126 (Springer, 2010).

Sémah, A.-M., Sémah, F., Djubiantono, T. & Brasseur, B. Landscape and hominid environment: Java (Indonesia) changes between the early Middle Pleistocene and the early Pleistocene. Quaternary ammonium salt. internationality. 223, 451–454 (2010).

Janssen, R. et al. The stable isotopes of enamel in the Holocene and Pleistocene fossil fauna reveal the glacier and interglacial paleoenvironment of ancient humans in Indonesia. Quaternary ammonium salt. science. Revision 144, 145-154 (2016).

Rizal, Y. et al. 117,000-108,000 years ago, Homo erectus last appeared in Ngandong, Java. Nature 577, 381–385 (2020).

CAS PubMed article Google Scholar 

Chen, F. etc. Mandibles of Denisovans from the Late Mid-Pleistocene from the Qinghai-Tibet Plateau. Nature 569, 409–412 (2019).

CAS PubMed article ADS Google Scholar 

Sutikna, T. etc. Revision of the Floresian stratigraphy and chronology of Liangbua, Indonesia. Nature 532, 366–369 (2016).

CAS PubMed article ADS Google Scholar 

Louys, J. & Roberts, P. Environmental drivers of the extinction of megafauna and humans in Southeast Asia. Nature 586, 402–406 (2020).

CAS PubMed article ADS Google Scholar 

De Vos, J. Rethinking the Pleistocene cave fauna of South China and its relationship with the Java fauna. Kur. Forsch. Research Institute Senkenberg 69, 259–266 (1984).

Schwartz, JH, Long, VT, Cuong, NL, Kha, LT & Tattersall, I. A variety of human fauna from the Late Mid-Pleistocene breccia cave in Tham Kuyen, Socialist Republic of Vietnam. Humanity. dad. Yes. Moose. Nat. history. 74, 1-11 (1994).

Schwartz, JH, Long, VT, Cuong, NL, Kha, LT & Tattersall, I. A review of the Pleistocene human fauna of the Socialist Republic of Vietnam. Humanity. dad. Yes. Moose. Nat. history. 76, 1-24 (1995).

Reyes-Centeno, H. Out of Africa and into Asia: Fossil and genetic evidence of modern origin and spread. Quaternary ammonium salt. internationality. 416, 249–262 (2016).

Bae, CJ, Douka, K. & Petraglia, MD On the origin of modern humans: an Asian perspective. Science 358, 9067 (2017).

Dennell, R., Martinón-Torres, M., Bermúdez de Castro, J.-M. & Xing, G. Population history of late Pleistocene China. Quaternary ammonium salt. internationality. 559, 4-13 (2020).

Westaway, KE et al. Early modern humans appeared on the island of Sumatra 73000-63000 years ago. Nature 548, 322–325 (2017).

CAS PubMed article ADS Google Scholar 

Bacon, A.-M. Wait. Late Pleistocene mammalian assemblages in Southeast Asia: New dating, mortality, and evolution of the predator-prey relationship in an environmental context. Paleogeography and palaeoclimate. Paleo ecology. 422, 101–127 (2015).

Bourgon, N. etc. Zinc isotopes in late Pleistocene fossil teeth from cave environments in Southeast Asia retain paleo diet information. The process of the national team Akkad. science. United States 117, 4675–4681 (2020).

CAS PubMed PubMed Central Article Google Scholar 

Bacon, A.-M. Wait. MIS6-5 Rhino megafauna during the transition period: Late Mid-Pleistocene Coc Muoi Assemblage, Langshan Province, Vietnam. Quaternary ammonium salt. science. Rev. 186, 123–141 (2018).

Bacon, A.-M. Wait. Nam Lot (MIS 5) and Duoi U'Oi (MIS 4) Southeast Asian Sites: Zoological Archaeology and Isotope Evidence. Paleogeography and palaeoclimate. Paleo ecology. 512, 132–144 (2018).

Suraprasit, K., Jongauttchariyakul, S., Yamee, C., Pothichaiya, C. & Bocherens, H. New fossil and isotopic evidence of the Pleistocene zoogeographic transition and the hypothetical savanna corridor of the Thai peninsula. Quaternary ammonium salt. science. Revision 221, 105861 (2019).

Sun, F. etc. Pleistocene mammalian paleo-ecology and paleo-climate changes in South China: Evidence from stable carbon and oxygen isotopes. Paleogeography and palaeoclimate. Paleo ecology. 524, 1-12 (2019).

Demeter, F. et al. Anatomical modern people in Southeast Asia (Laos) 46 ka. The process of the national team Akkad. science. United States 109, 14375–14380 (2012).

CAS PubMed PubMed Central Article ADS Google Scholar 

Shackelford, L. etc. Tam Pà Ling, Laos. Other evidence of early modern human morphological diversity in Southeast Asia. Quaternary ammonium salt. internationality. 466, 93–106 (2018).

Petraglia, MD, Breeze, PS & Groucutt, HS Blue Arabia: Check the colonization and diffusion model. In Geological setting, Palaeoenvironment and Archaeology of the Red Sea (eds Rasul, NMA & Stewart, ICF) 675–683 (Springer International Publishing, 2019).

Cappellini, E. etc. The early Pleistocene enamel proteome from Dmanisi addresses the Stephanorhinus phylogeny. Nature 574, 103–107 (2019).

CAS PubMed PubMed Central Article ADS Google Scholar 

Welker, F. etc. The enamel proteome indicates that the great ape is an early divergent horse. Nature 576, 262–265 (2019).

CAS PubMed PubMed Central Article ADS Google Scholar 

Welker, F. etc. Tooth proteome of human ancestors. Nature 580, 235–238 (2020).

CAS PubMed PubMed Central Article ADS Google Scholar 

Wang, W. etc. A sequence of mammalian fossils from caves in the Bubing Basin in southern China, including human teeth. J. Humph. evolution. 52, 370–379 (2007).

Rink, WJ, Wei, W., Bekken, D. & Jones, HL Chronology of Giant Panda-Saber-toothed Tiger Fauna and Great Apes in Guangxi Province, Southern China. Quaternary ammonium salt. Reservoir 69, 377–387 (2008).

Norton, CJ, Jin, C., Wang, Y. & Zhang, Y. Rethinking the "paleo-Arctic-Eastern biogeographic boundary" of Quaternary China. Paleoanthropology in Asia: From Africa to China and beyond (Editor-in-Chief Norton, CJ and Braun, DR) 81-100 (Vertebrate Paleontology and Paleoanthropology, 2010).

Turvey, ST, Tong, H., Stuart, AJ & Lister, AM Holocene survival of late Pleistocene megafauna in China: a critical review of evidence. Quaternary ammonium salt. science. Rev. 76, 156–166 (2013).

Ma, J. et al. Isotopic evidence of the foraging ecology of the Asian elephant (Elephas maximus) in southern China during the Late Pleistocene. Quaternary ammonium salt. internationality. 443, 160–167 (2017).

Owen-Smith, RN Megaherbivores. The impact of oversize on ecology (Cambridge University Press, 1988).

Louys, J. & Meijaard, E. The paleoecology of the large animal habitat in Southeast Asia during the Pleistocene and a review of environmental changes in the region. J. Biogeography 37, 1432–1449 (2010).

Graham, RW Diversity and community structure of the Late Pleistocene mammalian population in North America. Journal of Zoology. Fen. 170, 181–192 (1985).

Graham, RW The spatial response of mammals to environmental fluctuations in the late Quaternary. Science 272, 1601–1606 (1996).

CAS PubMed article ADS PubMed Central Google Scholar 

Price, GJ fossil bandicoot (Marsupial Bandicoot) and environmental changes during the Pleistocene in the Darling Hills, southern Queensland, Australia. J. System. Paleontology. 2, 347–356 (2004).

Stewart, The Progressive Effect of JR Individual Species' Responses to Quaternary Climate Change: An Analysis of the British Mammal Population. Quaternary ammonium salt. science. Rev. 27, 2499–2508 (2008).

Faith, JT, Rowan, J. & Du, A. Early humans evolved in dissimilar ecosystems. The process of the national team Akkad. science. United States 116, 21478–21483 (2019).

CAS PubMed PubMed Central Article ADS Google Scholar 

Zeitoun, V., Chinnawut, W., Debruyne, R., Frère, S. & Auetrakulvit, P. Sustainability review of mid-Pleistocene benchmark sites including giant panda-sabertooth fauna: Proboscis perspective . Quaternary ammonium salt. internationality. 416, 12-26 (2010).

Jablonski, D. & Sepkoski, JJ Jr. Paleontology, community ecology, and ecological model scales. Ecology 77, 1367–1378 (1996).

CAS PubMed article Google Scholar 

Graham, RW Quaternary mammal community: the correlation between individual responses and non-similar fauna. In Paleobiogeography: Generating New Insights In Coevolution of the Earth and its Biota (eds Lieberman, BS & Stigall, AL) 141–157 (Palaeontological Society Paper, 2005).

Stewart, the evolutionary result of JR individualism's response to climate change. J. Evolution. biology. 22, 2363–2375 (2009).

CAS PubMed article Google Scholar 

Hofreiter, M. & Stewart, J. Ecological changes, range fluctuations and population dynamics during the Pleistocene. Curry creatures. 19, R584–R594 (2009).

CAS PubMed article Google Scholar 

Tougard, C. and Montuire, S. Pleistocene environment reconstruction and mammalian evolution in Southeast Asia: Focus on the fossil fauna of Thailand. Quaternary ammonium salt. science. Rev. 25, 126–141 (2006).

Zeitoun, V. etc. The dating, stratigraphy, and burial of the Pleistocene site at Ban Fa Suai II (Northern Thailand): Contributions to the study of paleontological diversity in Southeast Asia. install. Palentor. 105, 275–285 (2019).

Williams, JW and Jackson, ST Novel climate, non-simulation communities and ecological surprises. front. Ecology. environment. 5, 475–482 (2007).

Bennett, KD and Provan, J. What do we mean by refuge? Quaternary ammonium salt. science. Rev. 27, 2449–2455 (2008).

Leonard, JA, Wayne, RK & Cooper, A. Population genetics of brown bears in the Ice Age. The process of the national team Akkad. science. United States 97, 1651–1654 (2000).

CAS PubMed PubMed Central Article ADS Google Scholar 

Leonard, JA etc. The extinction of megafauna and the disappearance of specialized wolf ectoderm. Curry creatures. 17, 1146–1150 (2007).

CAS PubMed article Google Scholar 

Barnes, I., Matheus, P., Shapiro, B., Jensen, D. & Cooper, A. The dynamics of the Pleistocene population extinction of the Bering brown bear. Science 295, 2267–2270 (2002).

CAS PubMed article ADS Google Scholar 

Hofreiter, M. etc. Before the last glaciation, European mammals lacked systematic geography. The process of the national team Akkad. science. United States 35, 12963–12968 (2004).

Shapiro, B. et al. The rise and fall of bison in the Bering steppe. Science 306, 1561–1565 (2004).

CAS PubMed article ADS Google Scholar 

Rohland, N. etc. Population history of existing and extinct hyenas. Mole. biology. evolution. 22, 2435–2443 (2005).

CAS PubMed article Google Scholar 

Gilbert, MTP, etc. Intraspecific phylogeny analysis of Siberian mammoths using the complete mitochondrial genome. The process of the national team Akkad. science. United States 105, 8327–8332 (2008).

CAS PubMed PubMed Central Article ADS Google Scholar 

Orlando, L. et al. Use ancient DNA to revise the recent evolutionary history of equines. The process of the national team Akkad. science. United States 106, 21754–21759 (2009).

CAS PubMed PubMed Central Article ADS Google Scholar 

Campos, PF etc. Ancient DNA analysis ruled out human beings as the driving force behind the population dynamics of the late Pleistocene musk ox (Ovibos moschatus). The process of the national team Akkad. science. United States 107, 5675–5680 (2010).

CAS PubMed PubMed Central Article ADS Google Scholar 

Campos, PF etc. Ancient DNA sequences indicate that the mitochondrial genetic diversity of the saiga antelope (Saiga tatarica) has been largely lost since the Pleistocene. Mole. Ecology. 19, 4863–4875 (2010).

CAS PubMed article Google Scholar 

Lorenzen, ED et al. Late Quaternary megafauna reacted specifically to climate and human species. Nature 479, 359–365 (2011).

CAS PubMed PubMed Central Article ADS Google Scholar 

Log, L. et al. Ancient DNA indicates that the origin of modern wolves can be traced back to the expansion from Bering Island in the late Pleistocene. Mole. Ecology. 29, 1596-1610 (2019).

Lord, E. etc. Population stability and adaptive genomic characteristics of woolly rhinos before extinction. Curry creatures. 30, 3871–3879 (2020).

CAS PubMed article Google Scholar 

Lister, AM The influence of the Quaternary Ice Age on mammalian evolution. Phil. Translated by R. Soc. London. B 359, 221–241 (2004).

Barnosky, AD The impact of climate change in the Quaternary on the species formation of mammals. J. Mammals. evolution. 12, 247–264 (2005).

Stewart, JR, Lister, AM, Barnes, I. & Dalén, L. Refugia revisited: Individual responses of species in space and time. Process R. Soc. B 277, 661–671 (2010).

Pushkina, D., Bocherens, H., Chaimanee, Y. & Jeager, J.-J. Late Middle Pleistocene snake cave diet and stable carbon isotope reconstruction of paleoenvironment in northeastern Thailand. Naturwissenschaften 97, 299–309 (2010).

CAS PubMed article ADS Google Scholar 

Suraprasit, K., Bocherens, H., Chaimanee, Y., Panha, S. & Jeager, J.-J. Inferred updates in northeastern Thailand from stable isotope analysis of Khok Sung herbivore enamel and land mammal cenogram Ecology and climate in the late epoch. Quaternary ammonium salt. science. Rev. 193, 24–42 (2018).

Suraprasit, K. et al. Long-term isotope evidence for the diet and habitat breadth of goats from the Pleistocene to Holocene in Thailand: impact on the extinction and conservation of the Himalayan goral. front. Ecology. evolution. 8, 1-16 (2020).

Bocherens, H. et al. Pleistocene South Asian mammalian diet and habitat flexibility: impact on the fate of the giant fossil ape Gigantopithecus. Quaternary ammonium salt. internationality. 434, 148–155 (2017).

Stacklyn, S. etc. Carbon and oxygen isotope evidence for the differentiation of diet, environment and niche of giant pandas and related mammals in South China during the Early Pleistocene. Paleogeography and palaeoclimate. Paleo ecology. 468, 351–361 (2017).

Ma, J., Wang, Y., Jin, C., Hu, Y. & Bocherens, H. Stable isotope (C, O) analysis reveals the ecological flexibility and differential survival of Pleistocene Eastern saber-toothed tigers and elephants in Southeast Asia . Quaternary ammonium salt. science. Rev. 212, 33-44 (2019).

Farquhar, GD, Ehleringer, JR & Hubick, KT Carbon isotope identification and photosynthesis. Anu. Pastor plant biology. 40, 503–537 (1989).

van der Merwe, New Jersey and Medina, E. Canopy effect, carbon isotope ratio and food web in the Amazon. J. Archaeol. science. 18, 249–259 (1991).

Zazzo, A. et al. Using stable isotope ratios of enamel carbonate to study the paleo-diet and paleo-environmental changes of herbivores in Pliocene Chad. Paleontology 26, 294–309 (2000).

Dansgaard, W. Stable isotopes in precipitation. Tellus 16, 436–468 (1964).

Longinelli, A. Oxygen isotopes in mammalian bone phosphate: a new tool for paleohydrology and paleoclimate research? Geochemistry. Cosmic chemistry. Journal 48, 385–390 (1984).

CAS Article ADS Google Scholar 

Luz, B., Kolodny, Y. and Horowitz, M. Oxygen isotopic fractionation between mammalian bone phosphate and environmental drinking water. Geochemistry. Cosmic chemistry. Journal 48, 1689–1693 (1984).

CAS Article ADS Google Scholar 

Fricke, HC, Clyde, WC & O'Neil, JR The intradental variation of mammalian enamel δ 18O (PO4) is used as a record of the seasonal variation of continental climate variables. Geochemistry. Cosmic chemistry. Journal 62, 1839–1850 (1998).

CAS Article ADS Google Scholar 

Fricke, HC, Clyde, WC, O'Neil, JR & Gingerich, PD Evidence of rapid climate change in North America during the recent Paleocene thermal maximum: Oxygen from biological phosphates in the Bighorn Basin (Wyoming) Isotope composition. Earth planet. science. Wright. 160, 193–208 (1998).

CAS Article ADS Google Scholar 

Kohn, MJ, Schoeninger, MJ & Valley, JW Herbivore. Oxygen isotopic composition of teeth: effects of diet and physiology. Geochemistry. Cosmic chemistry. Journal 60, 3889–3896 (1996).

CAS Article ADS Google Scholar 

Bryant, JD & Froelich, PN A model of oxygen isotope fractionation in the water of large mammals. Geochemistry. Cosmic chemistry. Journal 59, 4523–4537 (1995).

CAS Article ADS Google Scholar 

Kohn, MJ & Cerling, TE Stable isotopic composition of biological apatite. Priest minerals. Geochemistry. 48, 455–488 (2002).

Zheng, Z. & Lei, Z.-Q. Vegetation and climate change records of 400,000 years in the Leizhou Peninsula Volcanic Basin in southern China. Paleogeography and palaeoclimate. Paleo ecology. 145, 339–362 (1999).

Lee, S.-P. Wait. Pleistocene vegetation in Guangxi, southern China, based on palynological data from seven karst caves. Grana 59, 94–106 (2020).

Wang, Y. etc. The millennium and orbital scale changes of the East Asian monsoon over the past 224,000 years. Nature 451, 1090–1093 (2008).

CAS PubMed article ADS Google Scholar 

Chen, H. et al. The penultimate glacier monsoon record terminated by Huludong and two-phase glaciers. Geology 34, 217–220 (2006).

Article ADS CAS Google Scholar 

Kelly, MJ etc. The high-resolution features of the Asian monsoon between 146,000 and 99,000 BP in Donggak-dong, China, and the global relevance of events around Termination II. Paleogeography and palaeoclimate. Paleo ecology. 236, 20-38 (2006).

Milan, S. etc. Environmental conditions constitute the first evidence of modern humans in Tam Pà Ling, Laos: stable isotope records of carbonates from terrestrial gastropods. Paleogeography and palaeoclimate. Paleo ecology. 511, 352–363 (2018).

Bird, MI, Taylor, D. & Hunt, C. The paleoenvironment of Southeast Asian islands during the last ice age: the savanna corridor in Sundaland? Quaternary ammonium salt. science. Rev. 24, 228–242 (2005).

Marwick, B. & Gagan, MK Late Pleistocene monsoon variability in northwestern Thailand: The oxygen isotope sequence of the bivalve Margaritanopsis laosensis excavated from Mae Hong Son Province. Quaternary ammonium salt. science. Revision 30, 3088–3098 (2011).

Geist, V. On the relationship between social evolution and ecology of ungulates. Yes. zoo. 14, 205–220 (1974).

Bacon, A.-M. Wait. Test the savanna corridor hypothesis during MIS2: the Boh Dambang hyena site in southern Cambodia. Quaternary ammonium salt. internationality. 464, 417–439​​ (2018).

Cannon, CH, Robert, J., Morley, RJ & Bush, ABG Sundaland’s current refuge rainforest does not represent its biogeographic past and is highly susceptible to disturbance. The process of the national team Akkad. science. United States 106, 11188-11193 (2009).

CAS PubMed PubMed Central Article ADS Google Scholar 

Yuan, D. etc. The time, duration and transition of the last interglacial Asian monsoon. Science 304, 575–578 (2004).

CAS PubMed article ADS Google Scholar 

Habrin, J.-J. How old is the oldest Homo sapiens in the Far East? The process of the national team Akkad. science. United States 118, e2101173118 (2021).

CAS PubMed PubMed Central Article Google Scholar 

Boivin, N., Fuller, DQ, Dennell, R., Allaby, R. & Petraglia, MD The proliferation of humans in different environments in Asia during the Pleistocene. Quaternary ammonium salt. internationality. 300, 32–47 (2013).

Pereira, N. et al. People in the ancient rainforest: Late Pleistocene foragers at the Batadomba-Iena rock shelter in Sri Lanka. J. Humph. evolution. 61, 254–269 (2011).

Roberts, P., Boivin, N., Lee-Thorp, J., Petraglia, M. & Stock, J. Tropical forests and Homo. evolution. Anthropologist. 25, 306–317 (2016).

Roberts, P. & Petraglia, MD Pleistocene Rainforest: Obstacles or attractive environments for early human foragers? Archaeologists of the world. 47, 718–739 (2015).

Wedage, O. et al. About 45,000 years ago, Homo sapiens conducted a special rainforest hunting. Nat. Community. 10, 739 (2019).

PubMed PubMed Central article CAS Google Scholar 

Buck, G. et al. The "Human Revolution" in the Low Tropical Regions of Southeast Asia: The anatomical antiquity and behavior of modern humans in Niah Cave (Sarawak, Borneo). J. Humph. evolution. 52, 243–261 (2007).

Piper, PJ & Rabett, RJ Hunting in the rainforest: Evidence from the end of the Pleistocene in Lobang Hangus, Sarawak. internationality. J. Osteoarchaeol. 19, 551–565 (2009).

Mellars, P. Going East: New genetics and archeological perspectives on the colonization of modern humans in Eurasia. Science 313, 796–800 (2006).

CAS PubMed article ADS Google Scholar 

Posth, C. etc. The Pleistocene mitochondrial genome indicates a single major dispersion of non-Africans and a population turnover during the late glacial period in Europe. Curry creatures. 26, 827–833 (2016).

CAS PubMed article Google Scholar 

Roberts, P. & Stewart, BA defines the "generalist expert" niche of Pleistocene Homo sapiens. Nat. Humph. behavior. 2, 542–550 (2018).

Zachwieja, AJ, etc. Use the niche model in the Australian test case to understand human land preferences in the Late Pleistocene. Quaternary ammonium salt. internationality. 563, 13-28 (2020).

Shea, JJ Homo sapiens is like Homo sapiens: behavioral variability and "behavioral modernity" in Paleolithic archeology. Curry anthropologist. 52, 1-35 (2011).

Sun, X.-F. Wait. Ancient DNA and multi-method dating confirmed the late arrival of anatomically modern humans in southern China. The process of the national team Akkad. science. United States 118, e2019158118 (2021).

CAS PubMed PubMed Central Article Google Scholar 

Martinón-Torres, M. etc. On the misidentification and unreliable background of the new "Human Tooth" in Fuyan Cave (China). The process of the national team Akkad. science. United States 118, e2102961118 (2021).

PubMed article CAS Google Scholar 

Timmerman, A. & Friedrich, FT Late Pleistocene climatic drivers of early human migration. Nature 538, 92–95 (2016).

Article ADS CAS Google Scholar 

Kealy, S., Louys, J. & O'Connor, S. The least cost path model indicates the proliferation of northern humans from Sunda to Sahul. J. Humph. evolution. 125, 59–70 (2018).

PubMed PubMed Central article Google Scholar 

De Deckker, P. etc. Marine isotopes in Oceania Phase 4: The peak of intact glaciers 65,000 years ago: global connections and impact on human transmission. Quaternary ammonium salt. science. Rev. 204, 187–207 (2019).

Clarkson, C. et al. The human occupation of northern Australia 65,000 years ago. Nature 547, 306–310 (2017).

CAS PubMed article ADS Google Scholar 

O'Connell, JF, etc. When did Homo sapiens first arrive in Southeast Asia and Sahur? The process of the national team Akkad. science. United States 115, 8482-8490 (2018).

PubMed PubMed Central article CAS Google Scholar 

Brain, CK hunter and prey? Introduction to African Cave Burial Science (University of Chicago Press, 1981).

Lucchini, V., Meijaard, E., Diong, CH, Groves, CP & Randi, E. A new phylogenetic view of the Southeast Asian wild boar (Sus sp.) based on mtDNA sequence and morphometric data. J. Zoll. London. 266, 25-35 (2006).

Sponheimer, M. et al. Do "savanna" chimps consume C4 resources? J. Humph. evolution. 51, 128–133 (2006).

CAS PubMed article Google Scholar 

Cerling, TE etc. The diet changes of large herbivores from 4 to 1 Ma in Turkana Basin, Kenya. The process of the national team Akkad. science. United States 112, 11467–11472 (2015).

CAS PubMed PubMed Central Article ADS Google Scholar 

Tejada-Lara, JV, etc. Comparison of the isotope ecology of mammals in the western Amazon rainforest. The process of the national team Akkad. science. United States 117, 26263–26272 (2020).

Kohn, MJ The carbon isotope composition of terrestrial C3 plants is used as an indicator of (paleo)ecological and (paleo)climate. The process of the national team Akkad. science. United States 107, 19691–19695 (2010).

CAS PubMed PubMed Central Article ADS Google Scholar 

Regarding the original isotope data of Coc Muoi (Lang Son Provincial Museum, Vietnam) and Duoi U'Oi (Peace Provincial Museum, Vietnam), provided by Dr. Nguyen Gia Doi from the Institute of Archaeology, Vietnam Academy of Social Sciences Research permit was issued to AMB on November 17, 2017 certificate. The Heritage Department of the Vientiane Ministry of Information, Culture and Tourism of the Lao People's Democratic Republic issued the Tan Kengnan Fauna Isotope Analysis and Research Permit (Letter No. 495) to AMB on November 28, 2016. Funding came from the research laboratory BABEL (FRE 2029 and UMR 8045) CNRS/University of Paris, AMB, Department of Human Evolution, Max Planck Institute for Evolutionary Anthropology, Leipzig, Germany, NB (contract number 0117/037). FW was supported by the Marie Skłodowska Curie Personal Scholarship (n°795569) and the European Research Council (ERC) under the EU Horizon 2020 Research and Innovation Program (grant agreement n°948365). EC is supported by VILLUM FONDEN (n°17649). The author would also like to thank Clément Zanolli for helpful comments and Catherine Yvon for editing the English text.

UMR 8045 BABEL, CNRS, Université de Paris, Faculté de Chirurgie dentaire, 1 rue Maurice Arnoux, 92120, Montrouge, France

Department of Human Evolution, Max Planck Institute for Evolutionary Anthropology, Leipzig, Germany

Nicholas Bergen and Jean-​​Jacques Habrin

Applied and Analytical Paleontology, Institute of Earth Sciences, Johannes Gutenberg University, Mainz, Germany

Evolutionary Genomics Group, GLOBE Institute, University of Copenhagen, Copenhagen, Denmark

Frido Wilk and Enrico Capellini

UMR 7209 Archéozoologie, Archéobotanique: Sociétés, Pratiques, Environnements, Muséum National d'Histoire Naturelle, CNRS, Paris, France

Dennis Fiorello, Olivier Tombret and Elise Dufour

Department of Anthropology and Paleoenvironment, Institute of Archaeology, Hoan Kiem District, Hanoi, Vietnam

Nguyen Thi Mai Huong & Nguyen Anh Tuan

Heritage Department, Ministry of Information, Culture and Tourism of Laos, Vientiane

Thongsa Sayavonkhamdy, Viengkeo Souksavatdy and Phonephanh Sichanthongtip

Institut des Sciences de l'Évolution de Montpellier, Université de Montpellier, CNRS, IRD, EPHE, Montpellier, France

Ecole et Observatoire des Sciences de la Terre (EOST Géologie), Institut de Physique du Globe de Strasbourg (IPGS) (CNRS/UMR 7516), Institut de Géologie, Université de Strasbourg, Strasbourg, Fr​​ance

Philip Tinner & Quentin Bosch

UMR 7362 Laboratoire Image Ville et Environnement, Institut de Géologie, Strasbourg, Fr​​ance

Department of Earth and Environmental Sciences, Traps' MQ Luminescence Dating Facility, Macquarie University, Sydney, Australia

Geological Archaeology and Archaeology Research Group, Southern Cross University, Lismore, Australia

Institute of Vertebrate Paleontology and Paleontology (IVPP), Chinese Academy of Sciences, Beijing, China

Spitteurs Pan, Technical Cave Supervision and Exploration, La Chapelle, France

Eric Suzoni and Sebastian Frangul

Natural History Museum of La Rochelle, France

Department of Biomedicine, University of Minnesota School of Medicine, Duluth, Minnesota, USA

Department of Anthropology, University of Illinois at Urbana-Champaign

Lundbeck Foundation GeoGenetics Centre, GLOBE Institute, Copenhagen, Denmark

UMR 7206 Eco-Anthropologie, Muséum National d'Histoire Naturelle, CNRS, Paris, France

French Academy, Chairperson of Paleoanthropology, Paris, France

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

AMB, NB and ED were designed and studied; NB, ED, OT and DF were used for sample preparation; AMB was analyzed for morphometric data; NB and ED were analyzed for isotope data; FW and EC were analyzed for proteomics data; TS, VS , PS, AMB, NB, POA, PD, JLP, KW, RJB, QB, ES, SF, EPE, AZ, LS, FD, LS, JJH are part of the LAOS project; NTMH, NAT, AMB, POA, PD, JLP, KW, FD are part of the Vietnam project; AMB, NB, ED, JJH, and FW wrote this paper with the contributions of all other co-authors.

The author declares no competing interests.

Springer Nature remains neutral on the jurisdiction claims in the published maps and agency affiliates.

Open Access This article has been licensed under the Creative Commons Attribution 4.0 International License Agreement, which permits use, sharing, adaptation, distribution and reproduction in any media or format, as long as you appropriately indicate the original author and source, and provide a link to the Creative Commons license And indicate whether any changes have been made. The images or other third-party materials in this article are included in the article’s Creative Commons license, unless otherwise stated in the material’s credit line. If the article’s Creative Commons license does not include material, and your intended use is not permitted by laws and regulations or exceeds the permitted use, you need to obtain permission directly from the copyright owner. To view a copy of this license, please visit http://creativecommons.org/licenses/by/4.0/.

Bacon, AM., Bourgon, N., Welker, F. etc. A multi-agent approach to explore the arrival, environment and adaptation of Homo sapiens in Southeast Asia. Scientific Report 11, 21080 (2021). https://doi.org/10.1038/s41598-021-99931-4

DOI: https://doi.org/10.1038/s41598-021-99931-4

Anyone you share the following link with can read this content:

Sorry, there is currently no shareable link in this article.

Provided by Springer Nature SharedIt content sharing program

By submitting a comment, you agree to abide by our terms and community guidelines. If you find content that is abusive or does not comply with our terms or guidelines, please mark it as inappropriate.

Scientific Report (Sci Rep) ISSN 2045-2322 (online)